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I. INTRODUCTION 



The processes of electron capture are under intensive investigation of the experimentahsts 
and theoreticians during the last decades. Still there are some disagreements between the 
experimental and theoretical results. 

One of the observed processes is the process of Radiative Electron Capture (REC) ac- 
coinpanied by the emission of one photon. There are many experimental data available 
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l|, |2[ . REC is the dominant electron-capture channel in fast collisions of heavy ions with 
light target atoms. This process does not depend strongly on the interelectron interaction. 

Interelectron interaction becomes important in the process of the capture of two electrons. 
There are two different types of the processes with the capture of two electrons by the Highly 
Charged lons(HCl): Double Radiative Electron Capture (DREC) and Radiative Double 
Electron Capture (RDEC). DREC is a two-step process in which two uncorrelated electrons 
are captured in one collision and two photons are emitted. RDEC is a one-step process 
where the momentum and the energy of two correlated captured electrons are converted 
into energy and momentum of one emitted photon. The processes of double electron capture 
were investigated experimentally in js-G] and theoretically in j?-!^. 

In experiments UH RDEC was organized as a process where two free (or quasifree) target 
electrons are captured into a bound state of the projectile, e.g., into empty i^-shell of an ion, 
and the energy emerges with only one photon. This process can be treated as the inverse 
process to the double photoionization. But in RDEC unlike the photoionization the bare 
nuclei should be used. The RDEC is a convenient tool to investigate the electron-electron 
interaction in the processes of the ion-atom collisions. 

The first RDEC experiment 3| was performed with 11.4 MeV/u Ar^^^ ions hitting a 
carbon foil at the UNILAC in GSI in 1994. To obtain as much as possible high rate of 
double capture in one collision, solid carbon target was chosen. In this experiment the 
probability of the RDEC process is very small. Experimental data give only the upper limit 
of the value of the cross-section of the process: about 5.2 mb. 

Another RDEC experiment {s] was performed at the heavy ion storage ring (ESR) in 
GSI. Bare U^^"'" ions with the energy of 297 MeV/u have been used in collisions with gas of 
Ar atoms. From the data obtained in this experiment we can only conclude that the value 
of the cross-section is also very small: less than 10 mb. The upper limit of RDEC process 
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was estimated to be significantly lower than measured previously [4|. 

Recent experiment j6| with 38 MeV 0^+ ions shows that there is some discrepancy be- 
tween the experimental data and theoretical prediction. In this experiment ions of the 
oxygen 0®+ move trough the thin carbon foil. 



The RDEC process was investigated theoretically in |7H10l|. In these works the calcula- 
tions were performed within the nonrelativistic theory. 

The process of double electron capture is governed by the interelectron interaction. We 
investigate radiative double electron capture (RDEC) by a bare nucleus followed by emission 
of a single photon 

2e-(e) + X(^)+ ^ X'^^-^^+ilsls) + 7H . (1) 

The initial state is presented by two incident electrons 2e~ with the same energies e and a 
bare atomic nucleus X^^^^ with the charge Z. 

The final state is given by two-electron ion in the ground state X^^~'^^'^(lsls) and a single 
photon 'j{co). We concern only radiation of one photon, first, because in the experiment [g] 
one photon is registered, and second, since the radiation of one photon gives the major 
contribution to the value of cross-section of the process. 

We consider a scenario where the momenta of the incident electrons are the same for the 
both electrons. This scenario corresponds to the experimental situation . The results 

of our calculations are compared with available experimental and theoretical data. 

II. APPLICATION OF THE LINE-PROFILE APPROACH 

The double electron capture is a nonresonant process. However for its description we 
will apply the LPA first developed for the resonant processes. The LPA appeared to be a 
convenient tool for the description of the nonresonant processes as well. 

Within the framework of quantum electrodynamics (QED) atomic electrons are interact- 



ing with the quantized electromagnetic field and with the quantum vacuum HI]. Accord- 
ingly, the set of electrons (together with the atomic nucleus) is not a conservative system, 
and the concept of the energy for this system needs to be carefully treated. Within the 
LPA [12] the energy levels are associated with resonances in the natural line-profile for the 
process of resonant photon scattering. In order to keep the characteristics of the energy 
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levels independent on the particular features of the process of scattering, the resonance ap- 
proximation is employed. The resonance approximation consists of description of resonance 
area of the natural line-profile by the Lorentz contour which is characterized by two pa- 
rameters: position of the resonance and its width. The energy levels are connected with 
the corresponding resonances. The energy and width of an energy level are determined by 
position of the resonance and its width within the resonance approximation. 

The amplitude of the process of photon scattering (the initial and final states are assumed 
to be the ground state) can be presented as a matrix element of a special operator {12! . 
This operator can be constructed with employment of the QED perturbation theory. In 
general, this operator depends on the photon frequency (w) and can be considered as a 
complex symmetric matrix (in some basis set) or as a quadratic form reducible to a diagonal 
form. Within the resonance approximation the eigenvalues of this matrix determine the 
positions of resonances and their widths. The eigenvectors of this matrix are used for 
calculation of the transition probabilities. If we consider the probability of a particular 
transition between two energy levels, we need to calculate the amplitude of this process. 
The amplitude is derived as a matrix element of the photon emission (absorption) operator 
(also constructed with employment of the QED perturbation theory) calculated on the 
eigenvectors corresponding to the initial and final states. Application of the LPA to the 



evaluation of energies and transition probabilities is presented in detail in 12h14|. We note 



that the technique developed in 12l-ll4l| can be used only for the bound electrons. 

The aim of the present work is the evaluation of the cross-section for electron capture. 
This process can be considered as a transition 

/ ^ F, (2) 

where the initial state (/) in case of REC process corresponds to the two electrons: a 
bound Is-electron and an incident electron, i.e., continuum electron. The final state (F) 
is represented by two bound electrons in the j-j coupling scheme configuration. In case of 
RDEC the initial state corresponds to the two continuum electrons and the final state is 
the same. Since the initial state contains continuum electrons, the LPA can not be applied 
directly to these processes. However, we can introduce an auxiliary bound electron system 
which properties are explicitly related to the properties of original system. 

We can consider the highly charged ion being confined within a sphere of a large radius 
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R. Then, all the energy spectrum becomes discrete. If the radius of the sphere is large, the 
wave function of election confined within the sphere (spherical box) and the wave function 
of electron (with the same energy) not confined within the sphere almost coincide. Let the 
electrons have the energy e > mgC^. Eigenvectors of the Dirac equation for the point nucleus 
are well known [11]. The asymptotics (r — )■ oo) of the Dirac wave function for the electron 
in continuum reads 



1pejlm{r) 
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— cos(pr + 0o(r)) (4) 
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sin(pr + 0/(r)) , (5) 

Tcp 

where \Cg\ = |C/| = 1 and (l)g{r), 4>f{r) are the functions smoothly depending on r = |r|, 
= f /r. The energy (e) and momentum (p) of the electron are connected as = mgC^+p^, 
where me is the electron mass, c is the speed of light. The continuum electron function is 
normalized to the energy delta function. For the large radius R and coordinate r the electron 
wave function for the ion enclosed within a box is given by its asymptotic: Eqs. (jlj), ([5]). 
Accordingly, the difference between the nearest (discrete) values of the momentum (Ap) 
is defined by one half of the oscillation period of functions in Eqs. (HI), (E]) at the border 
{r = R): ApR = it. Then, the difference between the nearest values of the energy (e) is 

A. = PAp=PI. (6) 

e eix 

The equations in this Section should be understood in the asymptotic sense, i.e., the equa- 
tions are correct up to the terms disappearing when i? — )■ oo. 

The wave function of electron confined within the sphere of radius R can be written as 



(ill 



i^eM = -^Mr)eiR-\r\), (7) 



where ipe{f) is given by Eq. ([3]) 



is a normalization factor {{ipenir) is normalized to unity), 6{R — \r\) is the Heaviside step 
function. The normalization factor A^e^ is considered in detail in Appendix |X] (see Eq. (lASp ). 



Note, that the function ipe^lr) goes to zero at any r when i? — )■ oo as 



i^e,{r) ~ (9) 



At a fixed radius R the function ipe^{r) describes a discrete energy leveL Technically, we 
can consider a resonance corresponding to this level. So, the LPA can be applied to any 
energy level of the HCI confined within a sphere of a finite radius. 

Reasoning by analogy with the system not confined within a sphere, i.e., with a system 
which has a continuous part of energy spectrum [11], instead of consideration of a single 
energy level {eji), we have to consider all the energy levels within some interval 6e = [ei, 62]. 
The number of levels within this interval is proportional tol/R (see Eq. iQ ) . The integration 
over an interval [61,62] in the continuous spectrum is equivalent to the summation over all 
the states {n) with the energy (e„) from the interval [61,62] in the discrete spectrum (if the 
ion is enclosed within a sphere of radius R): 

where function F represents some physical property {e.g., cross-section). If the radius R 
goes to infinity, the number of discrete states in the energy interval [61,62] goes to infinity 
and the width of the energy interval containing only one state goes to zero. Accordingly, if 
6e Ae, we can write 

F{e) = i- 1 de'F{e') = i-F(n) = {N^,rF{n) , (11) 

Ae 

where 6„ = 6 is the only discrete state inside the energy interval Ae. Thus, the transforma- 
tion from the continuous to discrete spectrum results in the substitution of the continuous 
spectrum wave function ipe by the function ipe^ and in an additional factor 1/Ae = {Nej^)"^ to 
the function F (cross-section), where Ae is the distance between the nearest energy levels. 

We conclude that the LPA can be generalized to the case of continuum electrons in the 
initial or final states. We can introduce an artificial bound electron state e^ described by 
the wave function ipf,^. The energies and the angular quantum numbers of the continuum 
electron state ipe and the bound electron state ipea are equal. If there is one continuum 
electron in the initial or final states, the amplitude calculated with functions ijje and the 



amplitude calculated with functions ipen are related like 15] 



[/e = lim N^^Ue, . (12) 

R—>-oo 



If there are two continuum electron in the initial or final states, the amplitudes are related 
like 

Ue,e2 = lim Ne,^Ne^j,Ue,j,e2n , (13) 
it— ^oo 

where Ne-^ is the normalization constant for the corresponding electron given by Eq. ([7j). 
In this paper we will consider electrons with equal energies, accordingly, we can set N^f^ = 
Ne^j^ = Nf.^^. We note, that the limit i? — )■ oo is equivalent to limit A^^^ — t- oo. 

So, we can generalize the LPA for calculation of the amplitude of the process of the 
electron capture. We employ the artificial bound electron states cr defined by Eq. ((71) and 
apply the LPA for calculation of the transition amplitude {Ur), i.e., for the system where 
the continuum electrons are substituted by the bound electrons en. The amplitude of the 
electron capture is given by Eq. f[T^ (if there is one continuum electron) or by Eq. f[T^ (if 
there are two continuum electrons). The limit i? — )■ oo can be evaluated numerically. 



III. TWO-ELECTRON WAVE FUNCTIONS 



The incident electron can be characterized by momentum (p) and polarization or spin 
projection (/i), and described by wave function ipp^{r). The energy (e), momentum and 
electron mass (me) are connected as e = ^Jp^ + m^, where p = \p\. It is also convenient to 
introduce the electron wave vector i/p = p/\p\. The wave function of incident electron is 
normalized like 



/ 



dr^|J;,,,{r)^M = {2nf5{p'-p)S,,, (14) 

= ^^5(e'-e)5(cos0' -cos 0)5(0' -0)(5^.^, (15) 

where the set {p, 6, (p) represents the vector p in spherical coordinates. This normalization 
corresponds to one particle per unit volume. 

The wave function of the incident electron (ipp^) can be expanded in the complete set of 
wave functions (ipejim) with the certain energy (e), total angular momentum (j), parity (Z) 



and projection of the total angular momentum (m) 11 | 



j de y^^apf,^ejimtpejim{r) ■ (16) 

jlm. 
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Introducing the scalar product 



u,,Sm = (17) 



where VLjim{i^p) is the spherical spinor [16] and v^{vp) is the unit spinor function, the scalar 
products ap^^^'jim can be presented in the form 

f27r)3/2 



a. 



PIJ,,£jlm 



-5(£-e)e''^-'a,^^,,,™. (18) 



The phase ip^ji is determined by the field of the nucleus [ll|. The functions ipsjimif) are 
normalized like 

J dr^^,.,i,^,{r)^ejim{r) = 5{e' - e)Sj'j5in5m'm ■ (19) 

The wave function describing two incident electrons with the certain momenta and po- 
larizations can be written as 

^Piw,P2/.2(n,r2) = -^det{V^p,^,(ri),V^p,^,(r2)}. (20) 

We suppose that the initial state of the system is given by two incident electrons with 
the equal momenta (p) and the opposite polarizations (/xi = — /X2)- Accordingly, the wave 
function of the initial state is 

^''"Xr-i,r-2) = -^det{V'p^,=i/2(ri),^p^,=_i/2(r2)}. (21) 
IV. CROSS-SECTION 



The amplitude of the process of electron capture (?/«/) is defined via S-matrix [ll| 



Sif = {-2m)5iEf - Ei)Uif , (22) 

where Ei, Ef are the energies of the initial and final states of the system. Then, the transition 
probability is given by [17 1 



1 dk 

dw,f = 271— ^\U,f\'6{Ef - E,)^ , (23) 

where Ei, Ef are the initial and final energies of the whole system. The factor l/V"^ cor- 
responds to the densities of the incident electrons, V is the reaction volume. The wave 
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functions of the incident electrons are normalized by the condition: one particle per unit 
volume. The factor dk/ {27r)^ gives the number of photon states with certain polarization 
and momentum within an interval dk per unit volume: dk/{n^^{2TT)^), n''^^ is the photon 
density. The emitted photon is described by momentum (fc), frequency {u = \k\) and po- 
larization (A). Normalization of the photon wave function {A = {A^, A)), corresponding to 
one particle per unit volume, is 

/"rf^^{M)+(^)^(fc',A')(^) = (2nf—5{k-k')5x,x'. (24) 
J 2t<j 

Accordingly, the photon density {nJ'^) is set equal to unity. 



Cross-section is connected with the transition probability f l23|) as 11 1 



da,f = ^ , (25) 

where j is the current of the incident electrons. This current is defined as j = Tfv, where 
rf = 1/V and v = p/e are the density and velocity of the incident electrons, respectively, in 
the rest system of the nucleus. 

In the experiments js-sj the RDEC is considered as a process where a bare nucleus goes 
through target atoms and captures two electrons with emission of one photon. In our theo- 
retical model this process is considered in the rest frame of the bare nucleus. Accordingly, 
the incident electrons are located in the target atom. The reaction volume for one incident 
electron is 

Vt 

V=^, (26) 

where Vt is the reaction volume the target atom, Zt is the number of electrons in the 
target atom. We introduce the reaction volume V; within this volume the incident electron 
interacts with the nucleus. If the system is enclosed into sphere of a large radius R, then 
the reaction volume for the target atom (see Fig. [1]) is represented by cylinder which cross- 
section area is St = ttB^ {Rt is the target atom radius) and the length is equal to 2R: 
Vt = 2RSt. The reaction volume for one incident electron is 

V = 2RS, (27) 

where S is the area of the cross-section of the reaction volume for one incident elctron: 
S = St/Zt. The volume V can be expressed via the normahzation constant Nej^ (see Eq. 



V = ^ . 28 

e 

Here, we took into account that the incident electrons have the same energy, and therefore 
their normahzation constants {Ng^) are equaL 

In this work we calculate the total cross-section of the electron recombination, what 
means the integration over the directions of the emitted photon (i/^) and summation over 
the photon polarization (A). Then, we suppose that the incident electrons have the same 
momentum (Up), hence, we can also average over the electron momentum direction (vp). 

It is convenient to make a decomposition Eq. f |T6l) of the continuum electron wave function 
with certain momentum (p) and polarization (yu) over the electron wave functions with 
certain energy (e), total angular momentum (j), parity (/) and projection of the total angular 
momentum (m). The initial state (two incident electrons) is decomposed over the two- 
electron functions in the j-j coupling scheme. 



Accordingly, the cross-section can be written down as 



15| 



aif = lim . ,„ 



j dVkdUp\Ui^kXsf (29) 



where the photon frequency (u) is defined by the energy conservation law. The factor e/p 
in the square brackets comes from the current of the incident electrons. The factor l/Air 
represents the average over the direction of the momentum of the incident electrons (^'p); we 
suppose that the momenta of the incident electrons are equal. The factor N^^^, according to 
Eq. (|T3|) . shows that in the amplitude the one-electron wave functions (ipenjim.) are normalized 
to unity (see Eq. ([7])). The last factor in the square brackets is the contribution of the 
volume given by Eq. f l28|) . The first subindex {{) of the amplitude represents the initial state 
Eq. (1^ . The other sub- indices represent the final state: the sub- indices kX describe the 
emitted photon, the subindex s = {JsMgUsJ sjsins2j 82^2) corresponds to the two-electron 
configuration in the j-j coupling scheme 



^jAfniiiZin2i2/2(r'i,r2) = N C'jm("^i, 



7712) 



mim,2 



^ i'4'nijihmi i'^l)'4^n2j2h'm2 '4^n2j2hm2 (^l) V'nijiiimi 

(^2)), (30) 

where C-j ^{7711,1712) are the Clebsch-Gordan coefficients. The normalization constant is 
equal to l/\/2 for nonequivalent electrons and to 1/2 for equivalent electrons. 
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Expansion of the one-electron wave function with certain momentum and polarization 
over the wave functions with certain total angular momentum and parity (see Eqs. f|T6l) . 
(ITS]) ) and integration over the direction of the momentum of the incident electrons yield 



a. 



if 



lim 



R 



:2n 



,3\2- 



ep 



^ ^<^RhheR02l2eRhh(-RiAlA / d^^k UjMtRjilieiij2l2,k:\sUjMeRj3UeRj3lA,^l) 

J,MJl<j2,h<l2j3<jA,h<lA 

where coefficients AjM€iijiheRj2i2eRj3heRjAiA defined in Appendix [Bl The last factor in the 
square brackets in Eq. ( l3TI) comes from the expansion of the one-electron wave function over 
the wave functions with certain momentum and polarization (see Appendix [B|). The matrix 
element of the recombination amplitude is calculated with two two-electron wave functions: 
{JMejijlefij'l') for the initial state and s for the final state. 



V. CALCULATION OF THE AMPLITUDE 



Following the notations employed in 12|], we introduce the photon emission matrix ele- 
ments 

Aif^ = I drUr)l'Af>^\k)Ur) (32) 
and the one-photon exchange matrix elements 

(fi) = / rfridr2^„i(ri)^„2(r2)7f'7f /^,/,2(^],rl2)^dl(rl)V^d2(r2) . (33) 



The indices Ui, di designate one-electron Dirac states, Dirac matrices 7^' act on the one- 
electron functions tpdiifi), respectively. The photon wave function A^^''^^ is defined by Eq. 
(l2l|l . Function If^-^^^{Q,ri2) looks like 

I,.,2{^,n2) = (34) 

1^12 

"lJ.ltJ.2 JQri2 _|_ O I I e 



ri2 dxi^ 8x2^ ri2 

x(l-5^,o)(l-5^2o), (35) 



if Coulomb gauge is employed, or 



I,.,2{^,r,2) = (36) 
r 12 
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if Feynman gauge is employed. Tensor g^^^^^ is the metric tensor; Sf^-^^^^ is the Kronecker 
delta; ri2 = |ri — r2|. Repeated indices imply summation. 

The amplitudes presented in Eq. ( ETj) are defined by Feynman graphs depicted in Fig. |2J 
These amplitudes are proportional to the following expressions 

(■ \ ^ (■ \ ^ ^U2nluindid2 (<i'7\ 
^1 — Z^Sl,n — 2^ — ^ , W'J 

6 = Y.^2,n = Yl — — . (38) 

Here, C,i corresponds the contributions of the left {i = 1) and the right {i = 2) graphs in 
Fig. |2|, respectively. The indices di,d2 correspond to the continuum electrons of the initial 
state. The indices Ui,U2 correspond to the bound electrons of the final state. The index n 
may correspond to any state of the Dirac spectrum. 

If we consider an ion enclosed within a sphere of radius R, then for the states correspond- 
ing to the continuum we can write (see Eqs. ([1]), Q) 

where Ae is the distance between two closest energy levels. The asymptotics {R oo) of the 
matrix elements Aud and IuiU2did2 are investigated in Appendix |Al Here, we will investigate 
the behavior of the terms (i = 1, 2) with various values of the intermediate electron state 
n when R — i- oo. 

If n belongs to the discrete part of the spectrum, the terms C,i,n contain two wave functions 
of electrons from the continuous part of the spectrum: di and d2 states. Employing Eqs. 
flXTSjl . flXTTj) for and Eqs. (jM]), flMOj) for ^2,n we can write 

~ ^, 1 = 1,2. (40) 

The denominators in Eqs. ( 1371) . ( |38l) do not contain any smallness in this case. 
If n belongs to the continuous part, then in general case 

~ 1 = 1,2. (41) 

Formulas (1A19I) . (lAllI) and ( lAlOl) . ( ]A20I) are used for ,^i^„ and ^2,n, respectively. The denom- 
inators in Eqs. (1371) . (138!) are supposed not to contain any smallness. Consider now three 
special cases when Eq. (HTl) is violated. The first case is when the energy of the intermediate 
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electron state n is equal to the energy of the incident electron (e„ = eg): then, application 
of Eqs. fET9l) . flATT]) yields 

The denominator in Eq. f l37|) does not contain any smallness. For description of the two 
other cases it is convenient to introduce a continuum electron e with the energy 

eg = ee + ee - ei^ . (43) 

The second special case is when the energy of the intermediate state n is equal to the energy 
of the electron e (£„ = eg), then ,^1 „ — 00. Here, the energy of the intermediate two-electron 
state {sis + Sn) is equal to the energy of the initial two-electron state (eg + eg). The last case 
is when 5„ ^ eg. If Sn is the next state to eg (i.e., e„ = eg ± Ae), then 

~ ^ . (44) 



Here, formulas ( IA19|) . (]A11|) are used for A^^n and Iuindid2i respectively. The denominator 



in Eq. fl37j) is set to Ae given by Eq. ([6]). Consider contribution of the states contained 
in the interval (eg, eg + (5e], where 5t is a small finite value which does not depend on R. 
The number [K ^ 5e/Ae) of intervals Ae composing the interval 5e is proportional to R. 
Accordingly, we get 

£nG(eg,eg+5e] k=l 

Note, that the terms given by Eqs. (|lQ])-(j42]),(j44]),(j45]) vanish (faster than R'^/"^ or faster 
than N~^) in Eq. (!3T|) when N^^^ — )■ 00. Accordingly, they do not contribute to the limit 
in Eq. ( ISTl) . The nonvanishing terms come from the second special case where e„ = eg. In 
this case the denominator in Eq. (1371) is equal to zero and the standard perturbation theory 
is not applicable. However, the two-electron configuration {eji,eFi)j representing the initial 
state and the two-electron configuration {ls,6n)j, where J is the total angular momentum, 
can be considered as quasidegenerate ones. 

We have introduced the artificial electron state (e) by the condition: eg + ei^ = ee + e^. The 
subindex R at e^j indicates that the corresponding wave function is normalized to unity over 
the sphere of radius R. The configuration (Is, e^j) j has the same energy as the configuration 
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{^r,^r)j- The contribution of these configurations can be calculated within the framework 
of the LPA. The configurations (ls,e/{)j and {eR,eR)j are considered as quasidegenerate 
ones. 

I — I 

Within the LPA [12^ we compose the matrix V 



V = 1/^°) + AV 



(46) 



The matrix V"*^°^ is a diagonal matrix and it includes the one-electron Dirac energies cor- 
responding to a certain configuration. The matrix AV includes the one-photon exchange 
corrections as well as other corrections which can be omitted here. The matrix V can be 
written in a block form 



V 



Vn AVu 
AV21 V22 



(47) 



where the block Vu contains matrix elements constructed on the configurations mixing with 
the reference state (the initial or final state). The mixing configurations define the set g. 
The block V22 contains matrix elements calculated with all the other configurations, the 
blocks AV12, AV21 contains matrix elements constructed on one configuration from the set 
g and one not included in the set g. 

Consider the set g including only two configurations {ls,eR)j and {eR,eR)j given by 
the two-electron wave functions of the noninteracting electrons in the j-j coupling scheme 
(\&^?'' - ^ and '^f'^ s , respectively), where J is the total angular momentum. Then, the 
block matrix Vn is 2 x 2 matrix 



Vu 



(48) 



{Vn)u (AVn)i2 
{AV^^)2l {V^)22 

Note, that the matrix AV is composed by one-photon exchange matrix elements which in- 
clude continuum electron wave functions, therefore, they vanish when i? — )■ 00. Accordingly, 
(^11)11 = (^11)22 and (AVii)i2 = (AVii)2i — )■ 0, when i? — 00. The eigenvectors of the 
matrix Vu are 

V2 \ 1 



V2\l 



(49) 



The eigenvectors ($) of the operator V (represented by the matrix V) can be constructed 
with employment of the perturbation theory yj,]. In the zeroth order of the perturbation 
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theory, the eigenvectors ($(e_R,efl)j and ${is,efl)j) are combinations of the two-electron wave 
functions of the noninteracting electron in the ?'-?' coupling scheme - n and "^^^^ ^ ) 

, = (6,)2^!°) , +n(bi)i¥^K , (50) 

{eR,eR)j V {eR,eR)j 'V ^''^^ (l^.e^jj V 7 

, ^ = 7^(62)2^!°) , + (62)1^!? - , • (51) 

The factor 77 = ±1 defines which eigenvector corresponds to {eR,efi)j configuration. It can 
be determined by the sign of the mixing element (Al^ii)i2. 

The wave function "^fjj^en); proportional to R^^ (or to N^^), and therefore its con- 
tribution vanishes in Eq. (!3T|) when N^^^ — )■ 00. The wave function \E'|'j'^ -^^^ is proportional 
to R~^/'^ (or to A*"^"^) because it contains only one continuum electron. Accordingly, the 
contribution of the eigenvectors $(£^,6^),/, ^{is,eR)j is given by the contribution of the wave 
function Note that the admixture of the (ls,eij)j configuration to the {eR^e^jj 

configuration leads to a growth of the amplitude with factor A^^^^. This growth is compen- 
sated by the additional factor 1/V in the transition probability for the RDEC process (see 
Eqs. (|23|) and (|28|)). 

We can conclude that, if we take into account the interelectron interaction (the one- 
photon exchange corrections in all orders of the perturbation theory) between the Is, 
and cr electron states, the amplitudes of the processes 

e + e — )■ (Isls) + 7(0;) and Is + e — )■ (Isls) + 7(0;) 

are connected as 

U[e + e] = r]U[ls + e] . (52) 
This equation is vahd up to the terms disappearing when A^^^ — )■ 00. 



VI. NUMERICAL METHODS 

The Dirac spectrum is constructed in our work in terms of B-splines jisl . \\\ . The ion 
is placed into a spherical box with the radius = 70/ (aZ) (in the relativistic units), 
where Z is the nuclear charge and a is the fine-structure constant. The B-splines used in 
our calculations are of the order 8 and we employed the grid with 60 nonzero knots. The 
generated electron spectrum is discrete and finite. 
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The eigenvector and the corresponding eigenvalue {i.e., energy), which is the closest to 
the energy of the continuum electron (eg), is replaced by the wave function of the electron 
confined within the spherical box of radius R (Eq. ([7])), and by the energy eg, respectively. 
The electron states of the generated spectrum, which are close to the substituted electron 
state (e„), are designated as e„_i and e„+i. Extension of the number of the knots and the 
size of the box reduces the effect of the substitution of the continuum electron state (eg) by 
the e„ state. After some variations and tests with different conditions (the number of knots 
and the number of states when we sum over Dirac spectrum) we choose the number of knots 
and the size of the box to be large enough not to influence the accuracy of the computing. 

In the present paper we consider the RDEC to the ground state (Isls) and to the low lying 
single excited states (the KL-shell). The contribution of the i^L-shell states is determined 
by contribution of the (ls2s)o configuration. The contribution of the other states does not 
exceed 1% of the total cross-section. We also consider only the main channel of the RDEC: 
capture with emission of the electric photons with J = 1. 

The LPA [12I, [1^ is employed for calculation of the amplitude of the RDEC. In the 
framework of the LPA we fix a set of electron states and construct a set g of all possible two- 
electron configurations in the j-j coupling scheme built on this set of electron states. We 
introduce ns, np-electron states as the electron states in the B-spline approximation which 
energy is closest to the energy of the continuum electron state eg. The electrons included 
into the set g are Is, 2s, 2p, ns, np, {n±l)s, (n ± l)p-electron states. The (r;,±l)s, {n±l)p 
are electron states in the B-spline approximation next to the ns, nj9-electron states. 

In the numerical calculations we construct the matrix V in a. special way to include 



the contribution of the mixing configurations 12|, |l5|. The matrix V is calculated with 
application of the QED perturbation theory 

V{uj) = +AV = + V'^^\uj) + ... . (53) 

In physical sense the u is the frequency of the scattered photon and the matrix V depends 
on the value of the u. The matrix V^^^ is composed with one-electron Dirac energies cor- 
responding to a certain configuration. The matrix V^^\u) includes the first order QED 
corrections, such as self-energy (SE) and vacuum polarization (VP) corrections and one- 
photon exchange corrections. In our work the matrix V{ijj) contains only V^^^ and V^^\u), 



16 



the last one includes only one-photon exchange corrections: 



uiU2d\d2 



U2 



^^2 l)«l 



(54) 



(see Eq. ( |33|) ). Within the framework of the LPA the contribution of the one-photon ex- 
change correction is taken into account to all orders of the perturbation theory for the 
configurations from the set g. 

The amplitude of the transition from the initial state / to the final state F with emission 



of one photon with the frequency Uq can be written as [12 1 

Ui^F = (S(cjo))$„$, , 



(55) 



where S(a;o) is operator of emission of the photon, $/ and are the eigenvectors of the 
matrix V{ijj) corresponding to the I and F states, respectively. The operator E^uq) is 



evaluated with employment of the QED perturbation theory (see [12|, ll5|). In zero order 
approximation this operator coincides with the photon emission operator (A*^'^"'^"^*). In this 
work we consider only the one-photon exchange corrections to the operator S. According 



to 



15| . it reads 



H(o) + H«+eO(a^ 



The zero-order matrix element is 



^uiU2did2 



Op /l(fco,Ao)* e 



(56) 



(57) 



where An°,{2°^* are the matrix elements of the emission operator which includes the photon 
wave function Eq. ( l24ll . 

It is convenient to write the matrix ^ in a block form 

AV21 K?^ + AV22 



V 



Vu 


Vl2 




V21 


V22 





(58) 



where the block Vu is composed entirely with the states from the set g and the block V22 
does not contain states from the set g. The blocks V12 and V21 are constructed with one 
configuration from the set g and with another one not included in the set g. 
The matrix Vu can be diagonalized numerically (non-perturbatively) 



T/diag 



B'VuB 



(59) 
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where B is an orthogonal matrix, is the transposed matrix. Since in general is a 
complex symmetric matrix, the matrix 5 is a complex orthogonal matrix 



B'B = I. (60) 

where / is a unit matrix [lij = 6ij) of the appropriate dimension. 
The eigenvector of the matrix V can be written as [l^ . 

= E ^^.".^i? + E(^^^o.^. ^^;:'^Lo) ^r + • • • , (ei) 

where are the eigenvalues of the matrix Vu and E^^ are the sums of the Dirac energies. 
The functions "^^^^ are the two-electron functions in the j-j coupling scheme. The indices 
kg, Ig run over configurations from the set g, while the index k runs over configurations not 
included in the set g, i.e., over all the other two-electron configurations. The first term in 
the right-hand side of the expression Eq. fl6T]) can be considered as the zero order of the 
applied perturbation theory, the second term corresponds to the first order. 

The cross-section is given by Eq. ( l3T]) . where the amplitude U enters as its squared 
absolute value. Employing Eq. ( I6T1) . the amplitude U can be presented as f/ = t/^^-* -|- 
U^^^ + . . .. Accordingly, the squared absolute value of U can be written as |f/p = |f/*^°)p + 
2 Re f/*^^^} + |f/''^'*p. The last term corresponds to the second order corrections and, in 
principle, can be omitted. Still we prefer to keep it. We consider the contribution of this 
term as an estimate of magnitude of the higher order corrections. 

To compare our results with experiment we use in our calculation the model of two 
electrons which are moving along the same direction with equal momentum. 

The calculation was performed with different gauges for the exchange-photon (the 
Coulomb and Feynman gauges) and emitted photons (the transverse and nontransverse 



gauges 



12|). A small deviation of the gauge invariance took place. It is explained by the 
fact that the set of Feynman graphs that we take into account is not gauge invariant beyond 
the lowest QED PT order. The magnitude of the deviation is determined by the magnitude 
of the higher order corrections. 
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VII. RESULTS AND DISCUSSION 



In our work we calculate the cross-section for double electron capture by a bare nucleus 
followed by emission of the photon. We calculate the cross-section for three different ex- 
periments and present our results in Tables IIIIHVIII We consider the scenario when two 
electrons are going along one line and have the same value of the momentum p. In our 
model the incident electrons are considered as Dirac continuum electrons. According to 
Eqs. (123|) . (!25|) the cross-section of the RDEC process depends on the volume of reaction. In 
the RDEC experiments the captured electrons are initially located on an atom. We consider 
two approximation for the experiments: 1) we suppose that the electrons are distributed 
homogeneously in the atom ((jRdec,a^. 2) -^^g neglect all the electrons of the atom except 
the i^'-shell electrons and suppose that the electrons are distributed homogeneously within 
the sphere of the i^-shell radius of the atom (^o"Rdec,k-j 

In Tables Ht IIT] we present the radii and areas of the cross-sections of the target atoms (see 
Section HV|) used for the calculation of the RDEC cross-section. The radius i?^ is the radius 
of the target atom. The area of the cross-section of the reaction volume for one electron 
{S^) is calculated as = tt{R^)^ /Z^, where is the charge of the nucleus of the target 
atom. The area is employed in the case of the first approximation. The radius is the 
radius of the i^'-shell of the target atom. In the case of the second approximation the area 
of the cross-section of the reaction volume for one electron (S^) is set to = 7r(i?^)^/2. In 
the case of the second approximation (^o"I^dec,k-j ^^^q expressions for the RDEC probability 
and cross-section Eqs. (^^, (EI]) get an additional factor (R^/R^)^ (ratio between the 
volume of the i^-shell and the volume of the whole target atom). Accordingly, the values of 

^RDEC.A g^^^ ^RDECK ^iffg^g ^ factor 



First, we consider the case of the RDEC process when the ion of oxygen with the energy 
38 MeV is hitting a carbon foil. The density of the carbon foil is of the order of 10^^ 
particles/cm^. The results of our calculation for the RDEC process 0^^+C we present in 
the Tables llllt HVl IVl In Table [TTll we present the data for the RDEC process to the ground 
state. The first column presents the experimental value of the cross-section, the second gives 
the results of the nonrelativistic calculation The last two columns present our results. 




(62) 
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Data for the cross-section of the RDEC process to the (ls2s) state are given in Table IIV[ 
The RDEC to the (ls2s) determine actually the contribution of the RDEC to the all single 
excited states of the i^L-shell since the contribution of the states higher than (ls2s) are 
quite small. The total cross-section of the RDEC for oxygen (o-f'-i^Ec; _ (j^-^^^c _j_ ^rdec-j jg 
presented in Table El The experiment of the RDEC in collisions of O^"*" ions with carbon 
reported in |6,|. though the most detaUed descript.ou of the experiment is presented m 
[20|. We note that the the experimental data [20| for the separate contributions of the (Isls) 
and (ls2s) configurations were defined with the use of the calculations [t explains 

good agreement between the ratio of these contributions defined in 2o| and in [?, 9|. 

The RDEC experiment for argon was performed with 11.4 MeV/u Ar^^+ ions hitting a 
carbon foil at the UNILAC in GSI in 1994 [3]. The thickness of the target was 4-10 yUg/cm^. 
In this experiment the probability of the RDEC process is very small. Experimental data 
give only the upper limit of the value of the cross-section of the process: about 5.2 mb. The 
results of the calculations for RDEC process Ar^'^^+C are presented in the Table IVll 

The RDEC experiment for uranium was performed at the heavy ion storage ring (ESR) 
in GSI jsl. Bare U^^"^ ions with the energy of 297 MeV/u have been used in collisions with 
gas of Ar atoms. The density of the gaseous Ar-target was 5 x 10^^ particles/cm^. The 
data obtained in this experiment provide only the upper limit of the cross-section: less than 
10 mb. The upper limit of RDEC process was estimated to be significantly lower than 
measured previously ^. The results of the calculations for RDEC process U^^"^+Ar are 
presented at the Table IVnl 

Comparison of the results for the RDEC cross-section reveals disagreement between the 
experimental and theoretical data. The disagreement between the experimental data and our 
results can be explained by the model, which we employed for description of the target used 
in the experiments j^, 4, y]. This model is rather rough. In the experiments the captured 
electrons are initially the bound electrons of atoms of either carbon foil 3, ^ or argon gas 
^. In our model the incident electrons are considered as Dirac continuum electrons. In 
particular, we do not take into account the bound energy of the target electrons. We also 
suppose that the electron density in the target atoms is homogeneous. The disregard of 
the bound energy should exaggerates the results. The assumption that the electrons in the 
target atom are distributed homogeneously should also change the results, however, it is 
difficult to estimate its influence. 
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The disagreement between the theoretical results obtained in {q] and our results is not 
clear. In principle, the models employed for description of the target electrons are rather 
similar. The work jo] presents nonrelativistc calculation, while our one is fully relativistic. 
However the relativistic effects can not explain the present disagreement. We note that the 
volume of reaction (where the incident electrons interact with the bare nucleus) is defined 
different in these calculations. In our calculation the reaction volume is a cylinder (see Fig. 

n 

[1]), while in work [9| the reaction volume is a sphere of radius Rt, where Rt is the radius 
of the target atom in its rest frame. We also note that the ratios between the cross-sections 
of the RDEC to the ground state (Isls) and to the i^'L-shell (ls2s) (see Tables [ml HV]) 
obtained in these calculations are different. According to our calculation the total RDEC 
cross-section is determined by the RDEC to the ground state, while the results of 

m 

predict that the main contribution to the RDEC is given by the capture to the excited 
states. 

The model of quasifree electrons employed in this calculation is rather rough for descrip- 
tion of the experiments [sl, ^, where relatively light bare nuclear {Z = 8, 18) move trough 
carbon foil {Z = 6). Experiments where the target atoms are much lighter than the bare 
nucleus (e.g., experiment ^|) are preferable. This model would be also good for description 
of experiments with electron beams. The small relative velocity between the bare nucleus 



and electrons and the presence of the magnetic field could enlarge the cross-section 



2l|. 
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Appendix A: Asymptotics of matrix elements 

The asymptotics of the Dirac wave functions for electrons in the continuum is given by 
Eqs. (I3])-(I5]). These wave functions are normalized to the Dirac delta function with respect 
to the energy (Eq. (I19p ). In case of the normalization of the wave functions of the continuum 
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electrons to unity over a sphere of radius R the normahzation integral reads 

{N^.f = J dr^t,iJr)AMr)- (Al) 

r<R 

Employing Eqs. ©-(IS]) and performing integration over the angular variables (ur) we get 

R 

{N.J' = I dr{Mr)\' + \Mr)\') (A2) 



R 

cos^(pr + (pg{r)) H sin^(pr + ) (A3) 

np Tip J 



R 

dr [ \ cos(2pr + 20„(r)) H sin(2pr + 20^ (r)) . (A4) 

\Tip 27ip 2Tip J 



The first term of the integrand does not depend on R, its contribution is proportional to 
R. The last two terms contain the sine and cosine functions, the absolute value of their 
contribution does not exceed a value not dependent on R. So, we can write 

(N^.r = -i^ + OiR-'))- (A5) 

Accordingly, the wave function of the continuum electron normalized to the Dirac delta 
function and the one normalized to unity over the sphere of radius R are connected by Eqs. 

Consider the asymptotics {R — ?■ oo) of the one-photon exchange matrix elements 
{IuiU2did2) given by Eq. (1551) . For this purpose we can restrict ourselves by the Coulomb 
part of the function /^^^^(fi, ri2), which looks like 

I,.,Mn,) = (A6) 

r 12 

It is convenient to employ the decomposition 



^ I^^^fcK^^rJ, (A7) 



— 



where r< = min(ri,r2), r> = max(ri,r2), Pk{x) is the Legendre polynomial, Urii^r2 is the 
scalar product of vectors Uj.^ = Vi/ri {i = 1,2). For investigation of the asymptotics we can 
also retain only the term with A; = in the decomposition Eq. flA7p . So, we get 

4i«2did2 ^ J j ^^2^«i(^l)^i('^2)^^di(n)^d2(^2) ■ (A8) 

ri<R r2<R 
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We assume that all the electron wave functions are normalized to unity over the sphere 
of radius R. Accordingly, every continuum electron wave function has got a normalization 
factor l/A^e^ and it is proportional to I/VR (see Eq. ([9])). We will investigate the asymptotics 
{R — > oo) of the matrix element Eq. ( 1A8I) for various electron states n. 

At first we will suppose that the electron states ui, U2 correspond to bound electrons 
(e.g., Is-electron state), the electron state di (we will designate it as n) can be any state 
of the Dirac spectrum, the electron state d2 describes a continuum electron. If the electron 
state n corresponds to a bound electron, then there is only one continuum electron (^2) 
containing factor 1/\/R in the matrix element. Accordingly, we get 

If the electron state n corresponds to a continuum electron, then there are two continuum 



electron (n and ^2) containing factor l/y/R in the matrix element 

1 

Imu2nd2 ~ • (AlO) 

Now, we will suppose that the electron states di, ^2 correspond to continuum electrons 
with the same energies equal to e, the electron state ui corresponds to a bound electron 
(e.g., Is-electron state), the electron state U2 (we will designate it as n) can be any state of 
the Dirac spectrum. If the electron state n corresponds to a bound electron, then there are 
only two continuum electrons (di and ^2) containing factor l/\fR in the matrix element. 
Accordingly, we get 

Iuindid2 ~ (All) 

If the electron state n corresponds to a continuum electron (different from the electron ^2), 
then there are three continuum electrons {di, c?2 and n) in the matrix element. It yields 

Iuindrd2 ~ ^^372- (A12) 

We have to select a special case when the electron state n coincides with electron 1^2 • The 
matrix element Iurndid2 can be written as 

dri ^^i(n)^di(n) 



'iiindi d2 

ri<R 



— [ dr2^^{r2)i'd2{r2) + [ dr2'^^{r2)—^d2{r2) 
ri J J r2 

r2<r\ r\<.r2<R 
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(A13) 



If the wave function V'n is equal to iljd2, then the last term in the square brackets is 

/ c/r2^+(r2)-V^rf,(r2) = / c/r2^+ (r2)-^rf,(r2) (A14) 
J r2 J ^ r2 

r\<T2<R r\<T2<R 

R 

= N,^' I dr, - {\g,2{r2)? + \U{r,)\') , (A15) 
J r2 

ri 

where the asymptotics of the g^^ and functions are given by Eqs. (jl]), (|5]). Employing 
Eqs. dMJ-dAi]) we get 

rfr2V^+(r2)-V^,,(r2) = ^(log(/?) + 0(i?°)) , (A16) 
r2 n 

ri<r2<-R 

where the logarithmic term is given by the first term of the integrand in Eq. (lA4p . The first 
term in the square brackets in Eq. (1A13P is proportional to 1/R. Accordingly, for the case 
when the electron state n coincides with the electron d2 the matrix element Iuindid2 reads 

^uindid2 ~ ^3/2 • [-^^'J 

Consider also the asymptotics of the photon emission matrix elements Aud given by Eq. 
(|32|) . We assume that all the electron wave functions are normalized to unity over the sphere 
of radius R. If the both electrons {d and u) describe the bound electrons then the matrix 
element A^d does not depend on R 

Aud ~ . (A18) 

If one of the electrons {d or u) corresponds to the bound electrons and the other one to the 
continuum electron then 

Aud - (A19) 
If the both electrons describe the continuum electrons, then 

K 

Appendix B: Angular integration 



Integration over the direction of the incident electrons momentum in Eq. (129|) can be 
performed analytically. Only the two-electron wave function of the initial state depends on 
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the direction of the momentum. We consider the following integral (see Eqs. ( 129|) . (12T|) ) 
f I 1 



2 



The one-electron wave functions with certain momentum and polarization can be expanded 
in series over the wave functions with certain total angular momentum and parity. Employing 
Eq. f|T6|) we get 



><^£ljlhml(n)V'£2i2/2m2(^^2)V^*3J3^3„^3(r•3)V^*4i4^4m4(^4) (B3) 

The coefficients ap^p^ejim can be written as Eqs. (ITSl) . (ITTl) . The integration over energies 
yields 



AC "t'pMpiJi'imi'^ "fpMpiJ2i2m2^ "fpMpajaisms'^ ^'pAip3,i4/4m4 

Consider separately the following integral 

h = y"c/l/p(-l)-l/2+Mpi(_l)-l/2+Mp3 

MplMp3 

^'^I^pfj.pljlhmi<^i,p];^^j2hrn2^i^pfj.p3j3hrn3'^upjip,^j4hm4 ' (■^^) 

Employing Eq. ( ITTl) yields 

h = Y j djyp C^jJnl{'^ll^^^l)C^^l^2{^l2.^2)C^■J^^^ 

t^pl^lp3 miim;2m;3m;4/ll/^2/i3M4 

X^iimn (^p)^2m;2 i^p)^hm3 i^p)^l4mu i^^p) 
X (^_]^^1+Mpi+A'p3 

After integration over the direction of the momentum and summation over the proiections 

n 

|16l | we get 



4tt 

Kk 
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X(-l)'l+^2 + V2C_]^y3+i4 + l/2 

K 32 jl I I K J4 J3 




1/2 h h ] y 1/2 /3 h 

xC'j^i\mim2)C'^i\m^mi) . (B7) 



Here, we introduced the 3j-symbols, 6j-symbols [16| and 



n(jl,j2,...,Jn) = V2ji + lv/2j2+l... V2jn + 1. (B8) 



Accordingly, Eq. (1B4|) can be written as 

\3\ 2 



jlj2j3jAhl2hUm\m2mimA Kk 

K 32 ji 





1/2 h h 

This equation can be written in the form 



' Kk jl<j2,h<l2,j3<j4,h<U 



3\ 2 



A^12 '^Kki^l^2)det{lp,j^i^rm{ri)lpej2l2m2{r2)} 

J2 Cfk\^zm,) det{C,3,3^3(r3)C,.,,,„,(r4)} , (BIO) 



X 

m\m2 

where = 1/2 for equivalent electrons and iV = 1/a/2 for non-equivalent electrons; n 
{jlm). Here, we introduced coefficients 

^Kkejihejihmh^Uh = ^n(ji, j2, js, ^4, ^1, ^2, ^3, h) 




K 32 M ^ ^4 is 

1/2 /i /2 J \ 1/2 /3 /4 

X(— l)'2+J2 + l/2|^_-|^y4+J4 + l/2gjipejji^+jip,j2;2-j¥Jej3i3--jl^ej4;4 _ (Bll) 
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TABLE I: Radii (in pm) and areas of the cross-sections of the reaction volumes for one incident 
electron (in kilobarn) for the C target atom 







pK 




133 


919 


14 


31 



TABLE II: Radii (in pm) and areas of the cross-sections of the reaction volumes for one incident 
electron (in kilobarn) for the Ar target atom 











188 


617 


3.0 


1.4 



TABLE III: Cross-section (in barn) for RDEC process 0®++C, (tW^^s^ contribution. 



Experiment 


Theory 


-fi?S? [20] 


RDEC rqi 


RDEC.A . , 1 
^(isis) [this work] 


RDEC.K n, . ,1 
^(isis) [this work] 


3.2(1.9) 


0.15 


0.55 


0.019 



TABLE IV: Cross-section (in barn) for RDEC process 0^++C, contribution. 



Experiment 


Theory 




-.RDEC ryl 


RDEC,A r,, . 11 

^(1.2.) [this work] 


RDEC,K r,, . ,1 

^(1.2.) [this work] 


2.2(1.3) 


0.105 


0.05 


0.002 



TABLE V: Total cross-section (in barn) for RDEC process O'^^+C, (j^^^^c contribution. 



Experiment 


Theory 


^RDEC ^20] 


^RDEC 9J 


^RDECA [^i^ig ^^^^ 


^RDEC,K j^^-g 


5.5(3.2) 


0.26 


0.61 


0.021 



TABLE VI: Cross-section (in millibarn) for RDEC process Ari'^++C, af^^i^ contribution. 



Experiment 


Theory 


^RDEC rol 


^RDEC 
'^(Isls) 


[9] 


RDEC, A . , 1 
^(isis) [this work] 


RDEC,K . ,1 
^(1.1.) [this work] 


< 5.2 


3.2 


120 


4.3 
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TABLE VII: Cross-section (in millibarn) for RDEC process U^^'^+Ar, c^^-j^? contribution. 



Experiment 


Theory 


_RDEC 


_RDEC rnl 
"(Isls) 


RDEC,A ni . , 1 
^(isis) [this work] 


RDEC,K n, . ,1 
^(isis) [this work] 


< 10 


2.5 X 10-2 


1.73 


0.31 X 10-2 



FIG. 1: In the experiments the bare nucleus is moving through fixed target atoms. The present 
calculation is performed in the rest frame of the bare nucleus. The system is enclose into sphere 
of the radius R. The bare nucleus is fixed in the center of the sphere. The cylinder presents the 
reaction volume for the process of radiative double electron capture. The area of the cross-section 
of the cylinder (S't = ttR'^) is given by the radius of the target atom {Rt)- The volume of the 
cylinder is Vr = 2RSt- The reaction volume for one electron is y = Vt/^Ti where Zt is the 
number of electron in the target atom. 
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FIG. 2: The Feynman graphs representing the first order interelectron interaction corrections to 
the process of electron recombination. The internal wavy line denotes the exchange by the photon 
between two electrons. The indices di, d2 correspond to the initial one-electron states of a system; 
til, U2 correspond to the final states, the index n corresponds to the intermediate one-electron 
states. 
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